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APPLICATION  OF  RAY  THEORY  TO 
LOW  FREQUENCY  PROPAGATION 


INTRODUCTION 

It  is  often  said  that  ray  theory  is  not  applicable  to  low  frequency  propagation 
in  the  ocean.  The  purpose  of  this  report  is  to  demonstrate  that  this  is  not  the 
ease.  If  the  word  "ray"  is  allowed  a more  general  meaning  than  that  used  in  the 
classical  sense,  then  ray  tracing  is  indeed  a useful  means  of  modeling  low  fre- 
quency propagation. 

Early  ray  tracing  programs  were  primarily  concerned  with  integrating  the 
ray  tracing  equations  of  the  next  section  accurately  and  efficiently.  It  is  shown 
that  the  effect  of  sound-speed  representations  on  the  computed  value  of  propa- 
gation loss  is  not  as  important  as  is  currently  believed.  The  most  recent  addi- 
tion to  practical  ray  tracing  programs  is  the  asymptotic  treatment  of  caustics. 

In  the  case  of  a horizontally  stratified  ocean,  the  integral  representation 
may  be  expanded  into  a multipath  series,  each  term  of  which  corresponds  to  a 
particular  ray  type.  Upon  integrating,  one  obtains  the  acoustic  pressure  along 
the  ray.  It  is  important  to  note  that  this  multipath  expansion  is  exact.  The  ac- 
curacy of  the  final  result  depends  on  the  method  of  solving  the  depth  dependent 
wave  equation  and  evaluating  the  ray  type  integrals. 

For  low  frequency  propagation  in  nearly  horizontally  stratified  oceans,  the 
method  of  horizontal  rays  is  recommended.  Here,  the  pressure  is  expressed  as 
a summation  of  normal  modes  weighted  by  amplitudes  satisfying  horizontal  ray 
tracing  equations. 


RAY  TRACING  EQUATIONS 

Several  years  ago,  the  state  of  the  'irt  was  described  in  Officer's^  book  on 
sound  transmission.  Then,  ray  tracing  involved  approximating  the  solution  of 
the  reduced  wave  equation  for  the  acoustic  pressure  P 

v2  p + (-^-)2  P - 0 

with  the  Wentzel-Kramers-Brillouin-Jeffreys  (WKBJ)  form 
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P = a exp  (iwl)  . 
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The  travel  time  T and  amplitude  a satisfy  chc  eikonal  equation 


(VT)‘ 


and  the  transport  equation 
respectively. 


V-  az  VT  = 0 , 


The  eikonal  equation  may  be  solved  by  using  the  method  of  characteristics 
for  first-order  partial  differential  equations.  Characteristic  curves,  better 
known  as  rays,  are  orthogonal  to  surfaces  of  constant  time.  They  satisfy  the 
ray  tracing  equations 


d 

ds 


i n dx\  a i 

s \c  ds/  <3x  c 


/I  dy\  3 1 

\c  ds/  5y  c 


d / 1 dy 
ds 


d /I  dz\  d 1 

ds  \c  ds/  dz  e 


a_  _1 
ay  c 

a i 

az  e 


dT 

ds 


'dx\2  /dy\2  I dz\ 2 


ds^  + (dsj  \ds/ 


Once  the  rays  have  been  found,  the  divergence  theorem  applied  to  the  trans- 
port equation  produces  the  geometrical  spreading  law  for  the  pressure  amplitude, 


m 


{ 
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In  \6 


The  subscript  zero  refers  to  a reference  point  (usually  1 yd  away  from  a 
point  source),  and  bo  is  the  cross-sectional  area  of  an  infinitesimal  ray  tube. 
The  intensity  satisfies  the  conservation  of  energy  law  along  a ray  tube.  It  can 
be  shown  that  pressure,  on  the  other  hand,  satisfies  the  law  of  reciprocity. 

In  the  case  of  a horizontally  stratified  ocean,  that  is,  the  sound  speed  and 
ocean  boundaries  are  independent  of  both  horizontal  coordinates,  the  rays  re- 
main in  a vertical  plane  and  obey  Snell's  law, 

1 dr 

— — -•  constant  , 
c ds 


where 


i’2  = i y* 


SOUND  SPEED  REPRESENTATIONS 

When  ray  theory  was  first  implemented  on  digital  computers,  the  primary 
concern  wa”  to  integrate  the  ray  tracing  equations  accurately  and  efficiently. 
Pedersen^  motivated  the  design  of  many  ray  tracing  programs  by  demonstrating 
the  fact  that  discontinuities  in  the  sound  speed  gradient  cculd  introduce  anomalies 
in  the  computed  value  of  geometrical  spreading  loss. 

This  effect  is  illustrated  by  fitting  the  sound  speed  profiles  of  figure  1 with 
piecewise  linear,  1 piecewise  quadratic,  2 and  cubic  spline^  representations. 
Differences  are  more  readily  seen  in  the  sound  speed  gradients  shown  in  this 
figure.  The  corresponding  ray  diagrams,  figure  2,  show  that  caustics  due  to 
discontinuities  in  the  gradient  of  the  piecewise  linear  fit  disappear  when  smoother 
sound  speed  representations  are  employed.  Propagation  losses  computed  ac- 
cording to  classical  ray  theory  tend  to  accentuate  this  effect,  but  consider  what 
would  happen  if  a ray  theory  generalized  to  treat  caustics  correctly  were  used 
instead.  Figure  3 indicates  that  the  effect  of  different  sound  speed  representa- 
tions is  insignificant  providing  that  each  representation  accurately  describes  the 
input  data  to  be  fitted. 
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As  far  as  ray  diagrams  and  propagation  losses  computed  according  to  clas- 
sical ray  theory  are  concerned,  past  experience  indicates  that  cubic  splines 
produce  the  best  representations  for  analytic  type  sound  speeds.  However,  as 
the  input  data  become  more  irregular,  the  curve  fitting  procedure  becomes  more 
difficult  to  automate.  1 A second  disadvantage  of  cubic  splines  is  that  the  corre- 
sponding ray  tracing  equations  cannot:  be  integrated  in  closed  form,  which  is  a 
process  that  can  be  accomplished  with  piecewise  linear  and  quadratic  fits. 


Many  of  the  statements  made  above  are  also  true  when  the  sound  speed 
varies  with  one  or  more  horizontal  coordinates  as  well  as  depth.  If,  for  example, 
the  input  data  are  fitted  with  triangular  planes,  the  ray  tracing  equations  may 
be  integrated  in  closed  form,  but  anomalies  due  to  discontinuous  gradients  arc 
possible. 


ASYMPTOTIC  TREATMENT  OF  CAUSTICS 


In  the  last  fewycars,  significant  advances  in  practical  ray  tracing  techniques 
involved  the  treatment  of  caustics  rather  than  improvements  in  curve  fitting 
algorithms.  The  problem  may  be  illustrated  when  the  sound  speed  decreases 
inversely  as  the  square  root  of  depth,  as  shown  in  figure  1.  We  see  that  the 
ray  diagram,  figure  5,  forms  a well  defined  caustic. 
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Figure  4.  Sound  Speed  and  Gradient  Studied  by 
Pedersen  and  Gordon 
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Figure  5.  Ray  Diagram  for  a 1-kyd  Source  Depth 
Computed  by  Using  the  Sound  Speed  of  Figure  4 


Pedersen  and  Gordon5  compared  the  olassieal  solution  (solid  line)  with 
Brekhovskikh's5  modified  ray  theory  (broken  line)  in  figure  6.  Classical  ray 
theory  prediets  an  infinite  intensity  at  the  eaustie  at  3159  yd  and  an  infinite 
propagation  loss  in  the  shadow  zone  to  the  right  of  the  eaustie.  Pedersen  and 
Gordon  explain  that  the  abrupt  change  in  loss  at  3130  yd  oceurs  at  the  ray  that 
grazes  the  ocean  surfaee.  The  modified  ray  theory  did  not  apply  to  the  left  of 
this  ray. 

The  above  remark  points  out  the  difficulty  of  applying  modified  ray  theories 
to  the  simplest  of  eaustie  geometries.  Additional  effeets  due  to  the  ocean  bound- 
aries, eusped  eaustics,  ete. , prevent  the  theory  from  being  applicable  every- 
where. One  ean  program  as  many  speeial  eases  as  practical  considerations 
suggest,  but,  more  often,  one  uses  classical  and  modified  theories  outside 
their  domain  of  validity.  Sinee  eaustie  corrections  are  usually  obtained  by 
ineluding  additional  terms  of  a high  frequeney  expansion,  errors  increase  as 
the  frequeney  decreases. 
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Figure  7.  Propagation  Loss  for  a 0.  8-kyd  Receiver  Depth 
Computed  More  Accurately  Than  That  for  Classical  or 
Modified  Ray  Theories 


Consider  what  would  happen  if  a theory  generalized  to  treat  caustics  correctly 
were  used  instead.  The  result,  figure  7,  indicates  that  there  is  no  discontinuity 
in  propagation  loss  at  the  ray  that  grazes  the  ocean  surface  and  also  that  classical 
ray  theory  appears  to  be  more  accurate  to  the  right  of  the  grazing  ray  than  to 
the  left.  Consequently,  modified  ray  techniques  should  be  exercised  witheaution. 


Spofford  was  one  of  the  first  to  implement  modified  ray  theory  in  a practical 
computer  program.  The  procedure,  based  on  the  work  of  Ludwig,  7 assumes  that 
the  reduced  wave  equation  has  an  asymptotic  solution  of  the  form 

p = exp  (iwT)j  gAi  (w2/s  p)  + — — Ai ' (w2/3p)  j 

( iwl/3  f 


subject  to  the  orthogonality  condition 


VT  • V p - 0 , 
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where  Ai  is  the  Airy  function  of  the  first  kind,  and  T,  p , g,  and  h are  to  be 
found.  Upon  substituting  this  ansatz  into  the  reduced  wave  equation  and  compar- 
ing similar  powers  of  frequency,  one  obtains 

T - (T,  i T.)/2  , 

2/3  p3/2  - (T,  - T_)/2  , 

g p l^4  (a,  ^ a_)/2  , 

, -1/1 

h - p (a,  i a_)/2  . 


m 

l 


As  illustrated  in  figure  8,  subscripts  i and  - refer  to  the  two  rays  that  touch 
and  do  no*,  touch  the  caustic  before  reaching  the  field  point,  respectively.  There- 
fore, all  the  quantities  appearing  in  Ludwig's  representation  may  be  expressed 
in  terms  of  the  travel  times  Ti  and  amplitudes  a ( of  classical  ray  theory. 
Brekhovskikh's  solution  lacks  the  term  involving  the  Airy  function  derivative,  a 
term  that  is  important  away  from  the  caustic.  As  a result,  Ludwig's  solution 
has  a larger  domain  of  validity. 


Figure  8.  Classical  Rays  Used  to  Compute  a 
Uniform  Asymptotic  Solution  at  a Caustic 
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LOW  FREQUENCY  PROPAGATION  IN  HORIZONTALLY 
STRATIFIED  OCEANS 

Most  of  the  figures  discussed  before  were  produced  by  a computer  program 
designed  to  model  acoustic  propagation  in  a horizontally  stratified  ocean.  For 
mediums  such  as  this,  the  acoustic  pressure  due  to  a unit  point  harmonic  source 
situated  at  (0,  0,  zg)  has  the  integral  representation 


■f 


P (r,  z,  z ; o> ) - w I U\J  (to\  r)  G (z,  zs;X,  to)  dX  , 


where  the  Green's  function  G satisfies  the  depth  dependent  wave  equation 


£d2/dz2  < w-  | c 2(z)  - X2jj  G (z,  zs;X,to)  - -25  (z-zg) 


i 


and  suitable  boundary  conditions. 

The  method  of  solution  used  here,  that  is  multipath  expansion  of  the  integral 
representation,  is  quite  old,  dating  back  nearly  40  years  to  Van  dor  Pol  and 
Bremmer. iS  Following  Leibiger  and  Lee,  9 the  Green's  function  is  expressed 
in  terms  of  two  linearly  independent  solutions  F+  of  the  homogeneous  depth 
dependent  equation.  The  solutions  F>  are  normalized  so  that  their  Wronskian 
equals  -2u)i.  Upon  expanding  the  denominator  of  G intc  a geometric  series, 
the  double  summation 


P (r,  z,  zs; u 


oo  A 

E E 

v 0 n-1 


Pn)  (F’  Z' 


is  obtained.  If  z < zs,  one  sees  that 

P^'  (r,  z,  zs;  t o)=l  iwJ0  (uX  r)  F_(z;X,w)  F,  (z s;X,w) 

“o 

4r  (X’W)  7 hot  (X’W)  dX> 

where  7sur  and  “Y hot  are  boundary  reflection  coefficients.  Other  tc mis  of 
the  series  arc  similar,  each  integral  representing  a particular  ray  type.  The 
first  four  are  illustrated  in  figure  9.  It  is  important  to  note  that,  so  far,  the 
solution  is  exact.  The  validity  of  the  final  result  depends  on  the  method  of 
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Figure  9.  The  Principal  Ray  Types 


solving  the  depth  dependent  wave  equation  and  evaluating  the  ray  type  integrals. 
If  WKBJ  and  stationary  phase  techniques  are  used,  respectively,  the  classical 
ray  theoretic  solution  is  obtained.  MurphylO  replaced  the  WKBJ  technique  with 
a Weber  function  representation  in  order  to  treat  the  two-turning-point  problem. 
Brekhovskikh  replaced  the  method  of  stationary  phase  with  an  Airy  integral 
modification  in  order  to  investigate  caustics. 

At  present,  our  propagation  model  uses  the  following  modified  WKBJ  ex- 
pression to  solve  the  depth  dependent  wave  equation: 


F4  (z;\,w) 


1/2  1/6 


’ exp  j MwQ  ( 7 ;o,  zfc;X)  t i?r/-lj 


g(z;X)  (m  {«2/3P(z;X)}  iiAi{w2/3p(z;X)}), 


Q (zot  zt;x)  = [ {<-•  2(n-x2 } 


1/2 


df  , 


\ 2/3 

JP(z;X) 

j-  , and  g (z;  X ) - 

a* 

1/2 
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where 


zQ  is  a suitably  chosen  reference  point 

z ^ is  a turning  point 

Bi  is  the  Airy  function  of  the  second  kind. 

9 /o 

Whenever  w“/0p  is  a moderate  to  large  negative  number,  Ff  reduces  to  the 
usual  WKBJ  representation 

Fj_(z;X,w)  |e'“  (z)  - X2  j e\p|_iwQ  (z0,  z;X)|  . 

Since  this  modified  representation  is  inaccurate  in  the  vicinity  of  double  turning 
points,  I\  are  arbitrarily  set  to  zero  whenever  they  occur.  Hopefully,  this 
will  only  affect  a small  interval  of  integration  and  wall  not  introduce  significant 
errors  in  the  final  result.  Murphy's  technique  offers  an  alternative  procedure. 

The  method  of  evaluating  the  ray  type  integrals  is  based  on  the  following: 

1.  Segment  the  interval  of  integration  into  suitably  chosen  subintervals. 

2.  Use  stationary  phase  to  evaluate  subintegrals  whenever  possible. 

3.  Integrate  the  remaining  eases  numerically. 

It  was  originally  though  that  the  numerical  integration,  although  lengthy 
when  compared  with  stationary  phase,  would  be  invoked  so  infrequently  that  its 
contribution  to  the  total  computer  running  time  would  be  inconsequential.  So  far, 
this  has  not  been  the  ease.  Hopefully,  the  running  time  will  be  reduced  eventually 
w'hen  the  integration  routine  is  made  more  efficient. 

Since  it  is  customary  to  give  computer  programs  names  so  that  they  may 
be  distinguished  from  others  performing  similar  functions,  the  program  used 
herein  is  called  CONGRATS  V,  where  CONGRATS  is  an  acronym  for  Continuous 
Gradient  Ray  Tracing  System.  Actually,  the  completed  program  will  predict  the 
performance  of  active  and  passive  sonar  systems  and  is,  therefore,  more  than 
just  a propagation  program.  As  shown  in  figures  1 through  3,  CONGRATS  V has 
the  option  to  invoke  several  ray  tracing  procedures.  The  propagation  losses  were 
obtained  by  adding  the  multipath  contributions  coherently.  CONGRATS  V also 
produces  a plot  of  propagation  loss  using  power  addition,  in  wliieh  case  the 
phases  of  the  individual  contributors  are  neglected. 
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A COMPARISON  OF  PROPAGATION  MODELS 


At  present,  the  suite  of  the  art  of  propagation  modeling  for  stratified  oceans 
may  be  illustrated  by  two  figures  compiled  by  E.  Jensen  of  NUSC.  (See  figures 
10a  and  lla. ) Both  compare  FFP,  11  FACT,  12  RAYMODE  9,  1:3  and  NISSM  11 14 
computer  predictions  for  50-lIz  propagation  in  the  Pacific.  The  choice  of  pro- 
grams included  in  the  comparison  was  mainly  of  convenience,  since  each  is 
available  at  NUSC',  New  London,  and  all  but  FACT  were  designed  there. 


Briefly,  the  Fast  Field  Program  (FFP)  utilizes  Fast  Fourier  Transforms 
to  evaluate  the  integral  representation.  The  Fast  Asymptotic  Coherent  Trans- 
mission Model  (FACT)  is  a constant  gradient  ray  tracing  program  incorporating 
sophisticated  low  frequency  modifications.  RAY  MODE  9,  the  latest  version  of 
the  series,  uses  ray  theory  to  determine  which  intervals  dominate  the  integral 
representation,  but  uses  normal  modes  to  compute  the  acoustic  amplitude.  The 
Navy  Interim  Surface  Ship  Model  (NISSM)  II  is  a continuous  gradient  ray  tracing 
program  designed  to  predict  the  performance  of  active  sonar  systems.  Ail  but 
FFP  have  the  option  to  combine  multipath  contributions  incoherently  as  well  as 
coherently,  and  all  but  FFP  use  alternative  procedures  for  surface  duct  propa- 
gation. 


As  a result,  the  first  case  (figure  10a),  which  is  dominated  by  surface  duct 
propagation,  will  show  the  greater  variability.  FACT  is  an  order  of  magnitude 
faster  than  NISSM  II  and  RAYMODE  9,  while  FFP  is  a good  deal  slower. 


Upon  adding  CONGRATS  V to  the  comparison  (figure  10b)  and  invoking  the 
coherent  phase  option,  one  sees  good  agreement  with  FFP.  Had  the  incoherent 
phase  option  been  invoked  instead,  CONGRATS  V would  have  agreed  with  the 
others. 


The  second  case,  figure  lla  is  dominated  by  convergence  zone  propagation. 
The  agreement  is  better  than  before  although  running  times  continue  to  differby 
orders  of  magnitude. 


Upon  adding  CONGRATS  V to  this  comparison  (figure  lib),  one  obtains 
reasonable  agreement  with  FFP.  It  is  unusual  to  find  agreement  among  models 
that  are  based  upon  different  theories  and  written  by  different  programmers. 
Unfortunately,  comparisons  are  not  always  this  good.  Hopefully,  all  discrep- 
ancies will  soon  be  eliminated  or  at  least  accounted  for. 


Frequency  burface  Duct  Propagation 
(From  E.  Jensen,  of  NL'SC,  unput: 


CONGRATS  V Low  Frequency  Surface  Duct  Propagation  for  a Pacific  Sound  Speed  Profile 


SOURCE  DEPTH  500  ft  RAYMODE  9 (INCOHERENT) 8.7  sec 

RECEIVER  DEPTH  300  ft 

FREGUENCY  50  Hz  FACT  (INCOHERENT) 1.3  sec 


quency  Convergence  Zone  Propagation  for  a Pacific  Sound  Speed  Profile 
(From  E.  Jensen,  of  NUSC,  unpublished  data. ) 


SOURCE  DEPTH  500  ft  CONGRATS  2 T 


CONG  [UTS 


■ m ns 
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HORIZONTAL  RAY  THEORY  FOR  NEARLY  HORIZONTALLY 
STRATIFIED  OCEANS 


At  long  ranges  or  in  shallow  water,  the  effects  of  horizontal  variations  in 
the  sound-speed  or  bottom  characteristics  are  often  large  and  not  readily  modeled 
by  any  of  the  techniques  discussed  previously.  A nearly  horizontally  stratified 
ocean  is  one  in  which  the  horizontal  variation  is  slow. 


This  rather  vague  notion  is  quantized  by  introducing  a small  parameter  e 
and  assuming  that  the  properties  of  the  medium  depend  on  the  horizontal  coor- 
dinates X,  Y only  through  the  combinations^ 


x - e X,  y = e Y . 


This  being  so,  let  us  seek  solutions  of  the  redueed  wave  equation  in  the  form 


P (x,  y,  z;  e)«  exp  j 0(x,  y)/ie|  E A„  (x,  y,  z)  (ie)1'  . 

y=0 


Each  A^,  in  turn,  is  assumed  to  have  the  form 


oo 

A„  (x,  y,  z)  - E a[/(x,  y)  (x,  y,  z)  , 
k-0 


wi  ere  the  \pk  are  orthonormal  eigenfunctions  of  the  depth  dependent  wave  equa- 
ti  >n 


d^k  2 2 

— t-  K (x,  y,  z)^R=  XR 

dz 


subject  to  the  appropriate  boundary  conditions. 


Upon  substituting  our  ausatz  into  the  reduced  wave  equation  and  comparing 
similar  powers  of  ie,  one  finds  that  the  phase  function,  e , satisfies  the  hori- 
zontally dependent  eikonal  equation 


+ \ (x,  y)  , 

w/ 


where  Xp  is  one  of  the  eigenvalues,  \k,  computed  above. 

This  equation,  like  the  ordinary  eikonal  equation,  may  be  solved  by  using 
ray  tracing  techniques.  Note, however,  that  all  depth  dependence  is  missing. 
The  pressure  depends  on  depth  only  through  the  vertical  eigenfunctions.  It  may 
also  be  shown  that  the  leading  amplitude,  a(®),  satisfies  the  conservation  law 


WC. 

m 
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A computer  program  based  on.  tin*  consideration  above  was  written  to  proceed 
in  two  stages.  The  Li  rsl.  determines  the  eigenvalues  and  normalized  cigenl'iinc- 
Lions  at  each  point  ot  a rectangular  grid  in  the  horizontal  plane.  Then,  during 
the  second  part,  a set  of  horizontal  ray  tracing  equations  is  integrated  ior  each 
eigenvalue,  and  the  contributions  of  individual  modes  arc  combined  to  obtain 
the  total  field. 


As  in  ordinary  ray  programs,  only  the  leading  term  of  the  asymptotic  ex- 
pansion oi  each  mode  is  found.  The  expansion  I hen  reduces  to  that  derived  by 
Pierce  almost  10  years  ago.  1** 


The  program  predicted  propagation  loss  along  a 1300-iimi  track  extending 
northward  from  27  JIO'.N,  137  aO'W  to  32  :>()' N,  1:37  3'TW.  eleven  equidistant 
velocity-depth  profiles  obtained  from  the  measured  data  displayed  in  figure  12 
were  entered  into  the  computer  program.  Note  thal  the  SO  PAR  axis  rises  from 
a depth  of  71)3  m at  27  JIO'N  to  about  30  m at  32  UO'N.  Pack  of  relevant  data 
prevented  the  inclusion  of  any  dependence  of  sound  speed  or  bottom  depth  upon 
longitude. 


Figure  Pi  displays  propagation  losses  from  dynamite  charges  detonated 
300  ft  below  sea  level  along  the  track  to  a 2 300- ft  receiver  depth  situated  at 
27  MO'N.  The  top  graph  represents  observational  data,  while  the  middle  graph 
shows  computed  results.  The  two  are  superimposed  in  the  bottom  graph.  The 
figure  displays  an  interesting  feature*.  The  propagation  loss  elec  reuses  with  in- 
creasing range  beyond  12  N.  This  decrease  may  he  explainer!  by  Liu*  fiie-t  that 
the*  receiver  is  only  12  1 It  away  from  Ihe1  SO  PAR  axis,  where  the  signal  is 
strongly  affected  by  the  amplitude  of  the  lew  lowest  moeles,  as  shown  in  figure  1 I. 
As  the  source  ship  moved  north,  the  source  approache’d  IheSOPAR  axis  causing 
the  amplitude  of  these  modes  (figure  J3)  to  increase  to  sueii  an  extemL  that  even- 
tually the’  loss  due  to  horizontal  spreading  was  overcome  and  the’  total  propaga- 
tion loss  decreased. 


The  7 ,) , s rtf)— L'l.  receiver  depth  of  figure  Hi  is  well  below  the  turning  points 
of  the  first  few  modes,  and  so  the  signal  there  is  dominated  by  the  higher  modes. 
The  amplitudes  ot  these  modes  are  not  greatly  affected  when  the  source  approaches 
the  SO  PAR  axis;  therefore,  for  this  receiver,  cylindrical  spreading  dominates 
the  enti  re  track. 
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Propagation  Loss  versus  Range  for  a 2500-ft  Receiver 
a 500-ft  Source  Depth,  and  a 31-Hz  Frequenc}* 


FREQUENCY  31  Hz 
SOURCE  500  ft 
RECEIVER  10800  ff 


■i  i 1 i ■ 1—  * * * 

25°  N 30°  N 35°  N 40°  N 45°  N 50°  N 55°  N 60°  N 


Figure  16.  Propagation  Loss  versus  Range  for  a 10,  800-ft  Receiver  Depth, 
a 500-ft  Source  Depth,  and  a 31-Hz  Frequency 
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SUMMARY 


Contrary  to  popular  belief,  ray  tlmory  is  an  accurate  and  efficient  means 
of  investigating  low  frequency  acoustic  propagation  in  the  ocean.  Of  course,  in 
this  report  ray  theory  has  not  been  used  in  its  classical  sense. 


Several  illustrative  examples  proved  that  it  is  possible  to  design  a single 
propagation  model  that  agrees  with  analytic  solutions  and  measured  data,  as 
well  as  other  computer  programs,  i'lds  effort  is  more  difficult  than  one  may 
realize,  for  once  a computer  program  is  timed  to  the  actual  environmental  con- 
ditions of  a real  -_oan,  it  may  lx*  impossible  to  input  data  for  which  analytic 
solutions  are  known.  The  apparently  simple  task  of  comparing  programs  is  in 
reality  even  more  difficult  than  a comparison  with  analytic  solutions,  first, 
one  must  have  access  to  the  programs  being  compared.  Second,  the  programs 
must  treat  the  input  data  similarly,  finally,  the  programs  must  treat  the  out- 
put data  similarly.  For  example,  how  is  one  to  compare  coherent  nhasc  propa- 
gation loss  predictions  with  those  of  a random  phase  program? 


Although  all  the  computer  models  discussed  above  have  been  designed  within 
the  last  few  years,  the  theories  upon  which  they  are  based  arc  much  older. 
Therefore,  it  is  felt  that  improved  computing  facilities  rather  than  improved 
acoustic  theories  have  been  responsible  for  improved  prediction  capabilities. 


The  future  of  ray  theory  may  prove  quite  different.  Application  to  unstrat- 
ified media,  random  media,  etc.  is  the  next  logical  step,  but  these  theories  need 
more  development  before  they  can  be  implemented  into  practical  prediction 
models. 
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